One or more bound planets per Milky Way star 
from microlensing observations 
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Most known extrasolar planets (exoplanets) have been 
discovered using the radial velocity*"^ or transit^ meth- 
ods. Both are biased towards planets that are relatively 
close to their parent stars, and studies find that around 
17-30% (refs 4, 5) of solar-like stars host a planet. Gravita- 
tional microlensing* on the other hand, probes planets 
that are further away from their stars. Recently, a popula- 
tion of planets that are unbound or very far from their 
stars was discovered by microlensing^". These planets 
are at least as numerous as the stars in the Milky Way*". 
Here we report a statistical analysis of microlensing data 
(gathered in 2002-07) that reveals the fraction of bound 
planets 0.5-10 AU (Sun-Earth distance) from their stars. 
We find that 17^*% of stars host Jupiter-mass planets (0.3- 
10 Mj, where M]- 318 and is Earth's mass). Cool 
Neptunes (10-30 M^) and super-Earths (5-10 M^) are even 
more common: their respective abundances per star are 
52+^^% and 62+^^%. We conclude that stars are orbited by 
planets as a rule, rather than the exception. 

Gravitational microlensing is very rare: fewer than one 
star per million undergoes a microlensing effect at any 
time. Until now, the planet-search strategy^ has been 
mainly split into two levels. First, wide-field survey cam- 
paigns such as the Optical Gravitational Lensing Experi- 
ment (OGLE; ref. 11) and Microlensing Observations in 
Astrophysics (MOA; ref. 12) cover millions of stars every 
clear night to identify and alert the community to newly 
discovered stellar microlensing events as early as possible. 
Then, follow-up collaborations such as the Probing Lens- 
ing Anomalies Network (PLANET; ref. 13) and the Mi- 
crolensing Follow-Up Network (mFUN; refs 14, 15) moni- 
tor selected candidates at a very high rate to search for very 
short-lived light curve anomalies, using global networks of 
telescopes. 

To ease the detection-efficiency calculation, the observ- 
ing strategy should remain homogeneous for the time span 
considered in the analysis. As detailed in the Supplemen- 
tary Information, this condition is fulfilled for microlensing 
events identified by OGLE and followed up by PLANET 
in the six-year time span 2002-07. Although a number of 
microlensing planets were detected by the various collab- 
orations between 2002 and 2007 (Fig. 1), only a subset of 
them are consistent with the PLANET 2002-07 strategy. 
This leaves us with three compatible detections: OGLE 
2005-BLG-071Lb (refs 16, 17) a Jupiter-like planet of mass 
M ^ 3.8 Ms and semi-major axis a ^ 3.6 AU; OGLE 2007- 
BLG-349Lb (ref. 18), a Neptune-like planet (M ^ 0.2 Mj, 
a ^ 3 AU); and the super-Earth planet OGLE 2005-BLG- 
390Lb (refs 19, 20; M ^ 5.5 M^, a ^ 2.6 AU). 

To compute the detection efficiency for the 2002-07 
PLANET seasons, we selected a catalogue of unperturbed 
(that is, single-lens-like) microlensing events using a stan- 
dard procedure^', as explained in the Supplementary In- 
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Figure l: Survey-sensitivity diagram. Blue contours, expected number 
of detections from our survey if all lens stars have exactly one planet with 
orbit size a and mass M. Red points, all microlensing planet detections in 
the time span 2002-07, with error bars (s.d.) reported from the literature. 
White points, planets consistent with PLANET observing strategy. Red 
letters, planets of our Solar System, marked for comparison: E, Earth; J, 
Jupiter; S, Saturn; U, Uranus; N, Neptune. This diagram shows that the 
sensitivity of our survey extends roughly from 0.5 AU to 10 AU for plan- 
etary orbits, and from 5 Me to 10 Mj. The majority of all detected planets 
have masses below that of Saturn, although the sensitivity of the survey is 
much lower for such planets than for more massive, Jupiter-like planets. 
Low-mass planets are thus found to be much more common than giant 
planets. 

formation. For each light curve, we defined the planet- 
detection efficiency e(\ogd,\ogq) as the probability that a 
detectable planet signal would arise if the lens star had one 
companion planet, with mass ratio q and projected orbital 
separation d (in Einstein-ring radius units; ref. 22). The 
efficiency was then transformed^^ to e (log a, log M). The 
survey sensitivity S (log a, log M) was obtained by summing 
the detection efficiencies over all individual microlensing 
events. It provided the number of planets that our sur- 
vey would expect to detect if all lens stars had exactly one 
planet of mass M and semi-major axis a. 

We used 2004 as a representative season from the 
PLANET survey. Among the 98 events monitored, 43 met 
our quality-control criteria and were processed^"^. Most 
of the efficiency comes from the 26 most densely covered 
light curves, which provide a representative and reliable 
sub-sample of events. We then computed the survey sen- 
sitivity for the whole time span 2002-07 by weighting each 
observing season relative to 2004, according to the num- 
ber of events observed by PLANET for different ranges of 
peak magnification. This is described in the Supplemen- 
tary Information, and illustrated in Supplementary Fig. 2. 
The resulting planet sensitivity is plotted in blue in Fig. 1, 
where the labelled contours show the corresponding ex- 



2 



pected number of detections. The figure shows that the 
core sensitivity covers 0.5-10 AU for masses between those 
of Uranus /Neptune and Jupiter, and extends (with limited 
sensitivity) down to about 5 M®. As inherent to the mi- 
crolensing technique, our sample of event-host stars probes 
the natural mass distribution of stars in the Milky Way (K- 
M dwarfs), in the typical mass range of 0.14-1.0 Mq (see 
Supplementary Fig. 3). 

To derive the actual abundance of exoplanets from our 
survey, we proceeded as follows. Let the planetary mass 
function, /(log a, log M) = dA^/(dlogfl dlogM), where is 
the average number of planets per star. We then inte- 
grate the product /(log a, log M) 5 (log a, log M) over logo 
and log M. This gives E(f), the number of detections we 
can expect from our survey. For k (fractional) detections, 
the model then predicts a Poisson probability distribution 
P{k\E) - e^^E''/kl. A Bayesian analysis assuming an un- 
informative uniform prior P(log /) = 1 finally yields the 
probability distribution P(log f\k) that is used to constrain 
the planetary mass function. 

Although our derived planet-detection sensitivity ex- 
tends overalmost three orders of magnitude of planet 
masses (roughly 5 M® to 10 Mj), it covers fewer than 1.5 
orders of magnitude in orbit sizes (0.5-10 AU), thus pro- 
viding little information about the dependence of / on a. 
Within these limits, however, we find that the mass func- 
tion is approximately consistent with a flat distribution in 
log a (that is, / does not explicitly depend on a). The planet- 
detection sensitivity integrated over logo, or ^(logM), is 
displayed in Fig. 2b. The distribution probabilities of the 
mass for the three detections (computed according to the 
mass-error bars reported in the literature) are plotted in 
Fig. 2c (black curves), as is their sum (red curve). 

To study the dependence of / on mass, we assume that 
to the first order, / is well-approximated by a power- 
law model: /o (M/Mq)", where /o (the normalization fac- 
tor) and a (the slope of the power-law) are the parame- 
ters to be derived and Mq a fiducial mass (in practice, the 
pivot point of the mass function). Previous works on 
planet frequency have demonstrated that a power law pro- 
vides a fair description of the global behaviour of / with 
planetary mass. Apart from the constraint based on our 
PLANET data, we also made use in our analysis of the pre- 
vious constraints obtained by microlensrng: an estimate 
of the normalization'^ /o (0.36 + 0.15) and an estimate of 
the slope^^ (-0.68 + 0.2), displayed respectively as the blue 
point and the blue lines in Fig. 2. The new constraint pre- 
sented here therefore relies on 10 planet detections. We ob- 
tained 10-° ''2±"-22 (M/Mo)"°-^^*" '^ (red line in Fig. 2a) with a 
pivot point at Mq - 95 M^; that is, at Saturn's mass. The 
median of / and the 68% confidence interval around the 
median are marked by the dashed lines and the grey area. 

Hence, microlensing delivers a determination of the full 
planetary mass fimction of cool planets in the separation 
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Figure 2: Cool-planet mass function, a. The cool-planet mass function, 
/, for the orbital range 0.5-10 AU as derived by microlensing. Red solid 
line, best fit for this study, based on combining the results from PLANET 
2002-07 and previous microlensing estimates'^--^ for slope (blue line; er- 
ror, light-blue shaded area, s.d.) and normalization (blue point; error bars, 
s.d.). We find dW/Cd log a dlogM) = lQ-°-^^^°-^-(M/MsMr°-''^-°''' , where N 
is the average number of planets per star, a the semi-major axis and M 
the planet mass. The pivot point of the power-law mass function is at the 
mass of Saturn (A/sat = 95 Mj,). Grey shaded area, 68% confidence interval 
around the median (dash-dotted black line). For comparison, the con- 
straint from Doppler measurements-' (green line and point; error, green 
shaded area, s.d.) is also displayed. Differences can arise because the 
Doppler technique focuses mostly on solar-like stars, whereas microlens- 
ing a priori probes all types of host stars. Moreover, microlensing plan- 
ets are located further away from their stars and are cooler than Doppler 
planets. These two populations of planets may then follow a rather dif- 
ferent mass function, b, PLANET 2002-07 sensitivity, S: the expected 
number of detections if all stars had exactly one planet, regardless of its 
orbit, c, PLANET 2002-07 detections, k. Thin black curves, distribution 
probabilities of the mass for the three detections contained in the PLANET 
sample; red line, the sum of these distributions. 

range 0.5-10 AU. Our measurements confirm that low- 
mass planets are very common, and that the number of 
planets increases with decreasing planet mass, in agree- 
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ment with the predictions of the core-accretion theory of 
planet formation^^. The first microlensing study of the 
abundances of cool gas giants^' found that fewer than 
33% of M dwarfs have a Jupiter-like planet between 1.5^ 
AU, and even lower limits of 18% have been reported^''^". 
These limits are compatible with our measurement of S'^^'^o 
for masses ranging from Saturn to 10 times Jupiter, in the 
same orbit range. 

From our derived planetary mass function, we estimate 
that within 0.5-10 AU (that is, for a wider range of orbital 
separations than previous studies), on average 17^g% of 
stars host a 'Jupiter' (0.3-10 Mj), and 52+g% of stars host 
Neptune-like planets (10-30 M®). Taking the full range 
of planets that our survey can detect (0.5-10 AU, 5 M® to 
10 Mj), we find that on average every star has 1.6^^ ^^ plan- 
ets. This result is consistent with every star of the Milky 
Way hosting (on average) one planet or more in an orbital- 
distance range of 0.5-10 AU. Planets around stars in our 
Galaxy thus seem to be the rule rather than the exception. 
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Supplementary Information is linked to the online version 
of the paper at www.nature.com/ nature. 
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